Sensitive SQUID magnetometry for studying nano-magnetism 
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The superconducting quantum interference device (SQUID) magnetometer is one of the most sen- 
sitive experimental techniques to magnetically characterize samples with high sensitivity. Here we 
present a detailed discussion of possible artifacts and pitfalls characteristic for commercial SQUID 
magnetometers. This includes intrinsic artifacts which stem from the inherent design of the mag- 
netometer as well as potential issues due to the user. We provide some guidelines how to avoid 
and correct these, which is of particular importance when the proper magnetization of nano-scale 
objects shall be established in cases where its response is dwarfed by that of the substrate it comes 
with, a situation frequently found in the field of nano-magnetism. 

PACS numbers: 07.55Jg 
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I. INTRODUCTION 

Many developments in modern solid state and nano- 
scale physics make very sensitive magnetometry an es- 
sential tool to study the magnetic properties of samples 
like ultrathin films [l|, nanoparticles Q, nanodot arrays 
or dilute magnetic semiconductors (DMS) [3j which are 
all deposited on a substrate. Especially in the latter 
field recent theoretical predictions 
of anisotropic ferromagnetism [f| 

moments [1, Q or ferromagnetism in oxides like HTO2 
Q have triggered research efforts towards novel classes 
of material with magnetic order at room temperature 
and theories to explain the observed magnetic behavior 
[13, [ll|. Nevertheless, there are other studies which dis- 
cuss the influence of stainless-steel tweezers handling on 
the mag netic properties of HTO2 [H| and other diamag- 
nets [l3[, point out "possible pitfalls in search of magnetic 
order" for samples using sapphire substrates [14|, or dis- 
cuss potential contamination of the sample holder itself 
15]. Most of the exciting claims of unusual magnetic 
properties have the use of superconducting quantum in- 
terference device (SQUID) magnetometry in common. 
While some work on thin films is done using a home- 
built ultrahigh vacuum compatible SQUID magnetome- 
ter without external magnetic field [1], the majority of 
the experiments for solid state samples is performed us- 
ing the commercial SQUID magnetometer MPMS (XL) 
from Quantum Design [16]. Using the MPMS SQUIDs 
is widely spread mainly due to its high degree of user- 
friendly automation and reliability as well as the lack 
of commercial alternatives. Ironically, the (apparent) 
easiness of the MPMS's everyday operation, a virtually 
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effortless sample-to-" *.dat"-file transition, is largely re- 
sponsible for most of the erroneous reports. Despite the 
thick user-manual (three thick folders) the highly auto- 
mated measurement-routines switch off the indispensable 
critical discussion of the as-received data, an increasing 
bad practise in science, particularly among the recent 
breed of scientists who tend to accept computer-yielded 
results as given and practically error-free. This remark 
is particularly valid when the proper magnetization of 
nano-scale objects shall be established in cases where its 
response is dwarfed by that of the substrate it comes 
with. 

Here we present a detailed survey of possible artifacts 
found in measurements using high sensitivity magnetom- 
etry, like SQUID magnetometers, with the main empha- 
sis put on the MPMS-SQUIDs, as these have become now 
a kind of " laboratory standard" in nanoscale magnetism. 
We firstly discus issues stemming from the measurement 
technique in general and later proceed to artifacts which 
originate from the set-up of the measurement by the user. 
In parallel, we will provide an analysis of typical system- 
atic errors augmented with some recommendation how 
to avoid those potential research pit-falls. 



II. COMMERCIAL SQUID MAGNETOMETRY 

A. General remarks 

All measurements have been taken with MPMS XL 
(5 T) rigs at the University of Duisburg-Essen and Insti- 
tute of Physics of Polish Academy of Sciences in Warsaw. 
Nevertheless, the artifacts are visible in other data of low 
signal measurements recorded with various MPMS ma- 
chines, e.g., on Cr-doped InN [17[, GaAs single crystals 
[l8j or Gd ion implanted cubic GaN [n|, GaN doped 
with Fe [H[, GaN doped with Mn [H . Therefore they 
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FIG. 1. (Color online) An example of a typical "scan" col- 
lected by MPMS system during a single a measurement act 
using reciprocating sample option mode of measurement. The 
"Long Scaled Response" is the direct representation of V(z) 
without any compensation for system drifts applied. Its am- 
plitude corresponds to the strength of the measured moment. 
In this case the sample is GaN:Mn layer on 5x4x0.6 mm 3 sap- 
phire. This scan was acquired at H = 1 kOe and T = 15 K. 



can be taken as representative for the entire family of the 
MPMS instruments. For all measurements the RSO op- 
tion [2l[ (4 cm of sample movement and 0.5 Hz to 1 Hz 
of repetition frequency) of the MPMS was used which of- 
fers a higher sensitivity than the standard dc-transport. 
Whenever required we will use the unit emu, since the 
absolute size of the total moment m, irrespective of the 
sample volume, determines how close the actual magnetic 
signals are to the artefact-level of the instrument for low 
signal samples. The conversion from emu to the Si-unit 
A/m is straightforward: 1 emu= 10 3 A/m; the MPMS 
solely provides m in emu. 

Let us begin with a very basic statement that con- 
trary to their common name SQUID magnetometers do 
not measure directly any magnetic moment as e.g. a 
moving-iron instrument measures a flow of an electric 
current. Instead they merely read the current exited in 
the superconducting loop (called a flux transformer) by 
a change of a magnetic flux threading the loop, however 
even this is done in a technically elaborate way. The 
existence of the flux transformer allows to place the ac- 
tual SQUID sensor in a protected environment (usually 
in a niobium shield immersed in the liquid helium bath 
of the SQUID cryostat), sufficiently away from a sample 
chamber when the sample can be heated up to material 
limits (currently 800 K) and the magnetic field can be 
ramped up almost to the critical field of the supercon- 
ducting wire the flux transformer is made of. The de- 
tecting part of the flux transformer at the sample-end is 
wound outside a variable temperature insert of the mag- 
netometer, so it primarily transmit rather large changes 
of the magnetic flux caused by the change of the field 



in the external magnet (the source of experimental mag- 
netic field) and due to temperature and iJ-field variations 
of magnetization of all the components used to make the 
sample chamber (and actually everything else around). 
In order to minimize the spurious fluxes of the whole 
environment and further to separate the minute signal 
of the specimen from the former, two clever tricks are 
employed: (i) the sample-end of the flux transformer is 
shaped into a second order gradiometer, (ii) the sample is 
transported through the gradiometer sections producing 
a unique form of a time dependent current, </$(£), in the 
superconducting flux transformer. Now, as long as time 
variations of the other sources of the magnetic flux that 
couple to the flux transformer are sufficiently slower or 
weaker than the J$(t), a value of the sample magnetic 
moment m can be established through (most frequently) 
numerical comparison of expected and measured J$(t), 
with m being the proportionality factor scaling the for- 
mer into the latter. We remark here that obviously the 
sample travel across the detection coils takes place in a 
well controlled manner, so practically J$(i) is replaced 
by its dependence on the sample position in the sample 
chamber J$(z), converted to a "SQUID voltage" V(z) 
(which is usually called " scan" ) [22j . These types of scans 
are represented graphically by the SQUID software in its 
user interface (see Fig. [T] for the example). This char- 
acteristic shape of V{z) reflects the general form of the 
second order gradiometer, see the cartoon of its actual 
shape in Fig. [21 The expected gradiometer response func- 
tion T(z) [see Fig.^a)] is taken for a point dipole located 
on the z axis (the symmetry axis of the gradiometer). It 
should be noted that the assumption of an ideal point 
dipole limits the useful sample size to at most 5 mm 
along the sample travel to keep the systematic error in 
m below ^3%. 

Here we have to place a strong comment concerning 
the system sensitivity or resolution for absolute values of 
m. Both depend on the relevance of the established V(z) 
to the real experimental configuration and a fidelity of 
V{z) — ¥ m numerical conversion, so great care should be 
taken by the user on every step of the experimental proce- 
dure as well as later during analysis of the data provided 
by the magnetometer. This is a particularly important 
remark for everyone who investigates 'composite' speci- 
mens in which the signal in question constitute a small 
fraction of the total moment seen by the magnetometer. 
It can be either a (very) thin 'magnetic' layer deposited 
on the (typically) thick substrate, or a nano- or bio- ob- 
ject in a gelatin capsule or other container. It is a com- 
mon wisdom, as well as the manufacturer suggestion, to 
measure in such cases, at least once, the sample 'carrier' 
separately and subtract this information from the mea- 
surements performed on the whole assembly. 

In rare cases, the specimen is supported by a substrate, 
which has a pure diamagnetic response, i.e., field and 
temperature independent) diamagnetism, the substrate's 
contribution can be determined from the measured data 
of the composite sample itself and subsequently being 
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(i.e., field and temperature independent) diamagnetism, 
they can still be of the same order as the signal of the 
actual thin film or nano-object of interest. Typically, 
from the group of common substrate materials, only sap- 
phire or SiC fulfill the above criteria; however, great care 
is still in order to avoid magnetic contamination, which 
can originate at any stage in-between growth of the single 
crystal and mounting the sample for the magnetometric 
study (see below). 

Whichever way performed, this procedure usually leads 
to a high degree of "compensation", i. e. two large num- 
bers (with limited accuracy) of approximately the same 
size are subtracted to yield the desired small magnetic 
response of the actual sample in question. Thus, a high 
absolute precision of the data is a must. The key point 
in measuring weak moments is establishing of the abso- 
lute credibility of the data provided by the MSMS sys- 
tem. Before doing so, we remark here that throughout 
the paper we disregard the obvious influence of the ex- 
perimental noise because low as well as high frequency 
components can sufficiently be averaged out from V(z) 
by multiple scanning and position averaging of the data. 
Consequently, for the sake of simplicity, we treat the data 
as noise-free and of a sufficiently low standard deviation. 



B. Magnetic moment coupling 

1. Sample size effects 



FIG. 2. (Color online) Topmost cartoon sketches the second 
order gradiometer configuration. In the MPMS system the 
coil radius a — 0.97 cm and the section separation A = 1.519 
cm [2^] . (a) The gradiometer response to a point dipole travel- 
ing along z on the axis of symmetry. This is the actual shape 
of the reference T(z) scan that is used to compute sample 
moment m upon experimentally established scan V(z). (b) 
Calculated V(z) for the same point dipole as in (a) but 3 mm 
radially offset from z. (c) Same as in (b) but this time the 
moment of the dipole is located in xy plane. Bold black curves 
are calculated for the given situation indicated in each of the 
panels. Thin colored lines are least-squares regression fits of 
the reference T(z). The values of m indicate the relative in- 
crease (b) and reduction (c) of the magnetic moment reported 
by the system at these two experimental configurations. 

subtracted. However, this procedure is only adequate in 
cases where the substrate is devoid of paramagnetic im- 
purities (usually glass or MgO are full of these). Further, 
they should not have a small band-gap (like Si or GaAs) , 
since it gives rise to a temperature-dependent paramag- 
netic contribution stemming either from a van-Vleck type 
contribution due to the temperature dependence of the 
band-gap (l8j or from Pauli-like contributions stemming 
from temperature dependent carrier concentrations. Al- 
though these contributions are tiny for the substrate's 
overall magnetic response, after subtraction of the pure 



The point dipole-based approach to T(z) taken by 
Quantum Design (as well as other manufactures) has got 
profound implications for magnetometry as in the vast 
majority of measurements the specimens are up to 1 mm 
thick rectangulars of an area up to 5 x 5 mm 2 , which are 
not necessarily ideally symmetric and perfectly centered 
in the gradiometer. All these factors modify the shape 
of V(z) and by that influence the computed values of m. 
We start our analysis from size related effect and com- 
ment on the influence of the sample positioning in the 
next section. 

The problem of non-zero specimen dimensions and 
their influence on the strength of the coupling with pick- 
up coils (this affects both amplitude and the shape of 
the V{z)) has been recognized already very early [24]- 
[26j and the reader can find plenty of correction factors 
for various experimental configurations in these papers. 
However, these calculation are performed for objects of 
relatively high symmetry (spheres, squares, discs, etc.). 
For completeness, we list correction factors for cuboids of 
the most frequently used sample dimensions in Table |U 
In our calculations these bodies are approximated by a 
set of point dipoles evenly occupying the whole volume 
of the specimens and a superposition of fluxes created by 
each individual element to determine the shape of V(z). 
We also assume here that these bodies are uniformly mag- 
netized in the direction of the axis of the symmetry of the 
gradiometer (the z-axis). Since in such conditions V(z) 
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TABLE I. Correction factors for MPMS XL SQUID mag- 
netometers calculated for the most frequently met cuboid- 
shaped specimen sizes, a and b denote in-plane dimentions, 
w is the thickness. The correct value for the moment of an 
investigated specimen is obtained when the MPMS reported 
value for the "longitudinal moment" is divided by the relevant 
number from the table. 



Sample size 
(mm) 


Coefficients 
along the relevant side 




in- plane 


perpendicular 


a 


b 


w 


a 


b 


w 


3.0 


3.0 


0.0 


0.9938 


0.9938 


1.0113 


3.0 


4.0 


0.0 


0.9860 


0.9986 


1.0156 


4.0 


4.0 


0.0 


0.9891 


0.9891 


1.0219 


4.0 


5.0 


0.0 


0.9946 


0.9774 


1.0280 


5.0 


5.0 


0.0 


0.9831 


0.9831 


1.0341 


5.0 


6.0 


0.0 


0.9902 


0.9689 


1.0416 


6.0 


6.0 


0.0 


0.9758 


0.9758 


1.0491 


3.0 


3.0 


0.3 


0.9939 


0.9939 


1.0122 


3.0 


4.0 


0.3 


0.9985 


0.9847 


1.0170 


4.0 


4.0 


0.3 


0.9892 


0.9892 


1.0217 


4.0 


5.0 


0.3 


0.9951 


0.9775 


1.0278 


5.0 


5.0 


0.3 


0.9832 


0.9832 


1.0339 


5.0 


6.0 


0.3 


0.9902 


0.9690 


1.0414 


6.0 


6.0 


0.3 


0.9759 


0.9759 


1.0489 


3.0 


3.0 


0.5 


0.9940 


0.9940 


1.0120 


3.0 


4.0 


0.5 


0.9986 


0.9848 


1.0167 


4.0 


4.0 


0.5 


0.9893 


0.9893 


1.0214 


4.0 


5.0 


0.5 


0.9952 


0.9776 


1.0279 


5.0 


5.0 


0.5 


0.9833 


0.9833 


1.0336 


5.0 


6.0 


0.5 


0.9904 


0.9691 


1.0411 


6.0 


6.0 


0.5 


0.9760 


0.9760 


1.0486 



scans are expected to be symmetrical with respect to the 
center of the gradiometer (position zn = in Fig. [5]) con- 
trary to results presented in [26| both linear and iterative 
regression modes (the fitting fixed either at the center of 
the scan, or allowed to 'float' along z to chase the "peak 
position" , respectively) should give the same results (as- 
suming ideally centered sample at zq). The correction- 
factors which are smaller than unity for the "parallel" 
configuration originate from the vertical extent of these 
rather thin plates. Compared to an ideal point dipole 
centered within one of the windings of the gradiometer, 
a sample with a spacial extent along the gradiometer axis 
will induce a larger superconducting current of opposite 
sense in the other, more remote winding which leads ef- 
fectively to 'flattening' down the resulting V(z). A dif- 
ferent effect is responsible for correction-factors larger 
than unity for the "perpendicular" configuration, when 
the sample extends more in the xy plane than along the 
axis of the gradiometer. We know that the strength of the 
magnetic field produced by a current loop substantially 
increases when we move away from its center towards the 
rim. Thus, along the reciprocity principle, the coupling of 
a magnetic dipole with the sensing coil increases, result- 
ing in a 'taller' V{z). This case is exemplified in Fig.^b) 
where the test moment is placed 3 mm away from the axis 



of the gradiometer (r/a = 0.33). An apparent increase 
of to of about 6.5% will result from this particular case. 

Although the numbers listed in Table Q] fall to a range 
of small percents, and so they may be disregarded in 
rough test measurements, they are playing a decisive role 
in small moments determination, when a high degree of 
compensation is required. To illustrate this effect let us 
consider an orientation dependent studies of a magnet- 
ically weak and submicrometer thin layer deposited on 
a macroscopic diamagnetic substrate. DMS such as ni- 
trides or oxides doped with transition metals deposited 
on substrates are a good example of such samples. Let's 
further assume that we know the diamagnetic moment 
m^ a (H) of the substrate from a test measurement of a 
reference piece performed in only one orientation. Now, 
assuming 5 x 5 x 0.5 mm 3 dimensions of the specimen, 
the difference between parallel and perpendicular config- 
urations gives us about 5% difference in SQUID response 
(see Table U and [27j, Fig. 4). If the magnetic response as 
a function of orientation is then established by a simple 
subtraction of this reference m^ a (H), an iJ-dependent 
0.05 • m^i a {H) value will be artificially added (or sub- 
tracted) to the magnetic moments of the layer in question 
in the other configuration. Remarkably, since the satu- 
rated moments of these magnetically diluted layers are of 
the order of only few percent of the high field moments 
of their substrates, disregarding the effect of sample ge- 
ometry may dwarf their real magnetic response. 




Detection 
coil(s) 



FIG. 3. (Color online) Two in plane orientations giving the 
same difference in saturation moments for the edge contami- 
nant (the red star) as it would be between in-plane and per- 
pendicular orientations. Importantly, the radial offset r of the 
contaminant in (a) will be the same in the perpendicular case 
(not shown). 

Recently Garcia et al. [l3[ suggested that a point 
ferromagnetic contaminant located at an edge of a mil- 
limeters sized sample which changes its xy-position upon 
mounting in different orientations is responsible for large 
differences in saturation values, mg, seen in various 
orientation-dependent SQUID studies @, 0, IH, El- It 
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was argued that rotating the sample to measure in-plane 
and perpendicular orientations changes the position of 
the hypothetical contamination particle with respect to 
the axis of the gradiometer. This should yield as much as 
40% difference in mj values if that particle is located on 
an edge of a 3 x 3 x 0.5 mm 3 diamagnetic substrate, see 
Fig. 10 in 13]. Following that scenario the edge contam- 
ination was far from the symmetry axis of the gradiome- 
ter during perpendicular measurement, while mounted 
for in-plane measurement it found its place on the axis, 
as sketched in Figs.[5Jb) and (a), respectively. Although 
generally plausible, this effect can not produce a suffi- 
ciently large anisotropy in any MPMS system. The in- 
ner diameter of the sample chamber is about 9 mm, so 
the largest possible distance of a magnetic particle from 
the coils' axis is r = 4.5 mm (or r/a ~ 0.5) for which 
the increase of the coupling is only 17%. A 40% effect 
would require an r/a ~ 0.6, as calculated by the same 
authors, see Fig. 9, which is a location already outside 
the sample chamber. We finally note that for the sample 
studied in [H| r/a is 0.15 for contamination at center of 
an edge or 0.22 for location at the corner, which gives 
at most 2% difference for MPMS systems. Moreover, 
we find it highly fortunate to put this single edge con- 
tamination exactly at the magnetometer axis to see the 
maximum effect. In more realistic cases, as it is very 
hard to impose only a single contamination during han- 
dling, only smaller effects can be expected. The more 
evenly the contaminations are spread on the edges, the 
less difference in their mj will be seen between any two 
orientations. One should note, that the above effect is 
not only present in in-plane versus out-of-plane measure- 
ments. Also an in-plane rotation by 90° will create a sim- 
ilar effect, as it is exemplified in Fig. [3] but requires a less 
sophisticated sample mounting, i. e. additional contami- 
nation will be less likely. Interestingly, as it is elaborated 
in chapter IlIIl the experimental arrangement depicted in 
Fig. EJb) is very likely to produce a characteristic 'kink' 
in m(H) and a up-down jumps of the sample position 
reported by the MPMS at H where both substrate and 
contaminant signals compensate, see Fig. fTJl We finally 
note that for typical sample sizes in an MPMS system 
the effect of edge contamination can account for maxi- 
mum 6% difference in mg. Any larger anisotropy has to 
have a different origin. 

To conclude this part we remark that for accurate mag- 
netometry, the magnetometer response to the physical 
dimensions of the investigated bodies requires an accu- 
rate background compensation. In cases where the exact 
correction factors cannot be established numerically this 
requires reference measurements in basically all orienta- 
tions and geometries in question. 



2. Positioning of the sample within the gradiometer 

As long as the sample is mounted on the sample holder 
in a way that it is placed exactly on the axis of the gra- 



diometer, its exact position along the gradiometer axis is 
irrelevant for an accurate determination of m as it can 
always be centered either manually by the user or by a 
built-in automated routine (e.g., to compensate for the 
thermal expansion/contraction of the sample holder dur- 
ing temperature dependent studies). However, this is not 
the case when the sample is mounted off the gradiometer 
axis. As already mentioned above, the strength of the 
m O gradiometer coupling depends strongly on position 
of m with respect to center of the gradiometer coils, and 
so the correction factors listed in Table U are only valid if 
the specimen is ideally centered on the gradiometer axis. 
This centering is usually assured by the sample holder, 
and commonly used drinking straws (provided with each 
MPMS system by Quantum Design) can assure adequate 
centering of the material, as sketched in Fig [TUJa-b) . 

There are however at least three exceptions to this 
case. Firstly, drinking straws may not the best choice for 
accurate studies in nanomagnetism due to background 
signals, e.g., Ref. [l5| and other designs may not offer 
adequate centering, see chapter III Dl Secondly, even in a 
drinking straw the sample may not be placed centrosym- 
metrically in cases the sample is mounted as sketched 
in Fig HUT c). or a curved straw will displace the sample 
from the center. Finally, we note that in nanomagnetism 
the investigated material is usually deposited on one side 
only of a bulky carrier, and so it occupies a different vol- 
ume of the gradiometer than the substrate, and so a very 
precise absolute magnctomctry would require calculation 
of separate correction factors both for the substrate and 
for the layer. As the latter is very tedious, the best sug- 
gestion would be to investigate structures as similar as 
possible under very similar conditions, i.e. to use nearly 
identical and adequately tested sample holders and sub- 
strates. These allows to significantly reduce the number 
of control experiments with bare substrates and yet al- 
lows carrying out accurate magnetometry. 



3. Tilted moments (magnetic anisotropy) 

So far, we were considering only the case of vertically 
(along the z-axis) magnetized samples. A ferromagnetic 
anisotropy present in the investigated nanostructure (or 
less favorably existing in the whole specimen) , or a pres- 
ence of a transverse component of the external magnetic 
field can rotate m away from z. A transverse component 
of the magnetic field may appear in the superconducting 
magnet when it is near its remnant state, or when the sys- 
tem lacks a proper magnetic shielding so that the Earth's 
magnetic field, and/or a field produced by magnetized 
ferrous equipment or reinforced concrete bars close-by, 
can leak into to sample chamber. As the latter cases arc 
system- or environment-dependent, we leave them aside 
from more general considerations. On the other hand, 
the magnetic anisotropy itself is very frequently the ac- 
tual object of interest so this case will be considered in 
greater detail in the following. 
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FIG. 4. (Color online) Magnetization reversal process along the hard axis in a thin (d = 0), macroscopically large (5x5 mm 2 ), 
and uniformly magnetized (a single domain state assumed) layer exhibiting very strong uniaxial magnetic anisotropy as seen 
by a SQUID magnetometer equipped with 2nd order gradiometer. Experimental configuration is depicted in the main panel 
(a). The thick black line represents the z component while the grey line the radial component of m in this case. Both m and 
H are given in the relative units. Panels (b-f) depict the scans V(z) (green points) calculated at five characteristic values of 
H/ Ha- Then the ideal response function T(z) is fitted in both iterative (red dashed lines) and linear (solid blue lines) modes 
and the results are presented in the panel (a) as red diamonds and blue circles, respectively. 



In the axially symmetric gradiometer pick-up coils ev- 
ery magnetic moment can be decomposed into its verti- 
cal, m z , and radial, m r , components. Again, from the 
reciprocity principle, we know that m r will not couple 
to the SQUID as long as it is located anywhere on the 
symmetry axis of the pick-up coils (r = 0). On the other 
hand, although H r , the radial component of the magnetic 
field generated by a circular current loop, vanishes in the 
plane of this loop, m r (r ^ 0) will be always sensed by 
the SQUID's gradiometer due to both the vertical move- 
ment of the sample (the principle of operation) and a 
vertical separation of the constituent coils of the gra- 
diometer. This situation is exemplified in Fig. [U[c) indi- 
cating the general regularity that a presence of nonzero 
m r anywhere outside the gradiometer axis of symmetry 
produces an odd V(z). For a general case when both m z 
and m r are present at a certain volume inside the gra- 
diometer, V(z) acquires both even and odd components 
of a relative strength depending of the tilt angle of m 
and a degree of the magnetic mass distribution in the 
gradiometer coils plane (the xy plane). 

Let's illustrate this issue considering a magnetization 



reversal process along the hard axis in a thin [d = 0), 
macroscopically large (5x5 mm 2 ), and uniformly mag- 
netized (a single domain state assumed) layer exhibiting 
very strong uniaxial magnetic anisotropy. (Ga,Mn)As, 
a member of DMS family, is a good example of such a 
material. It can exhibit either perpendicular magnetic 
anisotropy (PMA) with the easy axis oriented along nor- 
mal to its plane, or in-plane magnetic anisotropy (IMA) 
with the normal being the hard axis, depending on the 
Mn content, hole concentration and temperature [30ll3lj]. 
As already argued before, stronger effects are expected 
when the sample is mounted in perpendicular configura- 
tion that is when the layer is parallel to the coils plane. 
We consider IMA case, so without extrnal magnetic field 
the magnetic moment of the layer is oriented in the sam- 
ple plane. Then m z — m z = 0, m r — |m r | = m, and the 
external field is needed to tilt and finally align m with the 
z axis when H becomes stronger than the anisotropy field 
Ha of this layer. For the uniaxial magnetic anisotropy 
m z = m for \H\ > Ha, and m z = m ■ H/Ha for - 
Ha < H < Ha , and a response linear in H is expected for 
w z . Fig. 2][a) depicts the experimental configuration and 
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FIG. 5. (Color online) Red: offset of the scan or the "sample 
position" with respect to the center of the gradiometer cal- 
culated in the iterative regression mode for the perpendicular 
case considered in Fig. [4] Brown: results for the in-plane case 
are given. The two horizontal dashed lines mark the border 
of the physical extent of a 5 x 5 mm 2 sample. 



shows the expected experimental response (bold black 
line). In panels (b)-(f) of the same figure green points 
depict calculated gradiometer responses in the presence 
of nonzero m z and m r at 5 different values H/Ha- Then 
we fit T(z) in both linear (solid blue lines) and iterative 
(dashed red lines) regression modes and plot the two ob- 
tained sets of m in the panel (a) as line-connected solid 
points of the same colors. 

Although at high fields the results of both methods 
coincide, one immediately notices a significant difference 
in the m values yielded by these two modes of fitting at 
low magnetic fields. Counterintuitive, the iterative mode 
which allows for reasonable good fits of V(z) in the whole 
range of magnetic fields yields progressively worse num- 
bers as the m r /m z ratio for this macroscopically large 
specimen increases. The red set of points reproduces 
neither the m r - nor the m z - dependence on H at weak 
fields. Actually, the magnitude of m reported by the it- 
erative regression mode for large m r /m z ratios depends 
strongly on the range of z$ values allowed for the fitting 
routine. At the same time, the linear mode, which does 
not fit the computed V(z) well at low magnetic fields, 
especially at H = 0, yields values of m which perfectly 
reproduce the expected values of m z (H) in whole range 
of H . In other words, the V(z) does not reflect solely the 
magnetic response of the specimen but also comprises 
artifacts stemming from the finite sample size and the 
gradiometric detection system. 

How can we know now that moment tilting is in force? 
If someone decides to use the iterative regression mode for 
m{H) loops at different temperatures then the signature 
of the maoment-tilting is included in the information on 
the fitted "best" sample position. This can be observed 
during the measurement or conveniently been read out 
afterwards from "Long Offset (cm)", the 7th column of 
the *.rso.dat file. Fig. [5] collects the "best" sample po- 
sitions of the iterative regression mode corresponding to 



the example presented in Fig. [4] The characteristic " di- 
vergence" around H = and a sudden change of sign 
indicate the presence of a substantial odd component in 
V(z). For comparison, Fig. [5] shows also the "best" cen- 
ter position for a similar sample investigated in-plane. 
Then largest m r /m z ratio will occur near H = when 
the system exhibits PMA or when a system with a strong 
in-plane uniaxial anisotropy is measured along its hard 
axis, or oriented at any direction at H near the coercive 
field when the magnetization reversal occurs through co- 
herent rotation of the whole m. Please note two features 
in this figure: (i) although the artificial off-centering is 
smaller for the parallel case, the divergence is equally 
strong in both cases; (ii) the apparent shift of the sample 
position during the m(H) loop exceeds the real sample 
dimensions. These features are observed despite the fact 
that the sample cannot change its position during the 
measurement at constant temperature. The only remedy 
in this case is to recalculate the raw data having the sam- 
ple position fixed for the whole m(H) at the real sample 
position, established, for example, from high field data. 
The only way to avoid recalculation is to set the system 
to work in the linear regression mode, but this require 
proper centering the sample before each field loop, so to 
compensate for the temperature expansion of the sample 
holder. In general, we may not retrieve 100% accurate 
values, but at least the odd component will not severely 
mar the outcome. Note, that the discussed here feature 
is very similar in appearance to the edge-contamination 
presented in Fig. IT27c). 

C. The magnetic field 

The ability to precisely control the values of the exter- 
nal magnetic field acting on the studied specimen is one 
of the most important features of any useful magnetome- 
ter. The MPMS XL system offers magnetic fields up to 
50 or 70 kOe which are generated in a high-homogeneity 
superconducting coil which is energized by a high current 
power supply controlled by the MPMS software. How- 
ever, there are some issues related with the generation of 
the magnetic field that influence the fidelity of the final 
outcome of the magnetometer. The most severe cases 
may drive the inadvertent scientist to false or even un- 
physical conclusions. These issues are related to both the 
magnet power supply and to the magnet itself, and will 
be discussed in that order. We also exemplify an effect 
of a substantial reduction of the MPMS accuracy seen 
at magnetic fields stronger than about 25 kOe. As we 
cannot unambiguously identify its origin, we will discuss 
it at the end of this section. 



1. Magnet power supply 

The contemporary power supplies offer a very high 
degree of linearity of their voltage-to-current (input-to- 
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output) conversion. The current for the magnet of the 
MPMS magnetometer is provided by a unipolar power 
supply (Kepco) . The very important point to be noticed 
by any user of the system is that the actual magnetic 
field is not measured directly by a field sensor. Instead, 
the field values are determined by the voltage-drop over a 
shunt resistor which is amplified by a factor of 21 before 
it is compared to the set voltage of the control board 
to select the desired magnetic field. The difference of 
the two voltages is fed into the Kepco as control voltage. 
This magnetic field control has to be well-calibrated on a 
regular base and procedures are available from Quantum 
Design both for the voltage amplifier and the zero-current 
offset of the Kepco. 

To cover the wide field range of the MPMS there are 
two shunt resistors with respective amplifiers between 
which the system is switching when is moving from high 
resolution to low resolution charging mode. The spe- 
cific field of switching is related to the field/current mag- 
net calibration factor (unless forced by the user-defined 
sequence), but always within region near few kOe. It 
is factory calibrated to better than +/- 0.5% but the 
exceptional precision of the SQUID and relative charg- 
ing accuracy of the Kepco enable the step to be clearly 
seen. In particular the switch is marked in Figs. |5] and 
[7] by faint dashed lines (this MPMS switches just above 
5 kOe). Ideally, the adequate calibration should remove 
the step, practically, one can only minimize it by ap- 
plying slightly different correction factors at both field 
ranges. The other, perhaps more practical, solution is 
to treat the slightly wrongly reported field values as an 
additional contribution to a remnant field in the magnet 
and correct both effects simultaneously, as described in 
the next section. 

It should be further noted that the magnetic field con- 
trol comprises a self-correcting routine which is intended 
to correct small discrepancies between the actual cur- 
rent in the magnet and the current set by the power 
supply which can stem from moderate flux relaxation. 
In cases of field losses of the magnet, malfunction of 
the persistent-switch, or severe flux creep in the mag- 
net, this self-correcting mechanism may lead to increas- 
ingly larger correcting factors resulting in very inaccurate 
nominal field values compared to the real ones. This can 
result in large residual hysteresis or even inverted arti- 
ficial hysteresis. To remove these issues, the 1822-unit 
(DSP board) of the MPMS has to be switched off for at 
least 10 sec (the message: '"power-on reset occurred"' 
has to be flagged). This procedure which resets the mag- 
netic field correction routine is also necessary before the 
zero-current offset of the Kepco is corrected. 



2. Superconducting coil 

As stated above, there is no absolute field sensor lo- 
cated anywhere in the magnetometer, so the magnetic 
field value reported to the user is based only on the cur- 



rent set by the power supply. Because a superconducting 
solenoid is used to generate the magnetic field, the user 
has to take into account the existence of pinned magnetic 
flux lines and flux movement within the magnet. In gen- 
eral, this makes it difficult to precisely determine the ac- 
tual magnetic field at the location of the sample. This is 
particularly true at weak magnetic fields. There are two 
phenomena that cause a discrepancy of the MPMS nomi- 
nal field values (which is laid down in the data files) from 
those really experienced by the sample. Namely they are 
superconducting magnet remanence and flux creep and 
escape [32] . In the following we will discuss the main im- 
plications of these effects for the studies of small-signal 
samples and suggest methods of mitigating these arti- 
facts. 

The remnant or trapped field in the superconducting 
magnet, H tr , results from the penetration of the material 
of the magnet by quantized magnetic flux. This flux gets 
pinned at defects in the material and when the magnet 
is discharged some of these flux lines remain pinned, giv- 
ing a weak magnetic field at the sample location. Impor- 
tantly, the spacial distribution of H tr is very complicated 
(see for example [33| , where apart of a characterization of 
this problem in MPMS systems, some methods of achiev- 
ing the lowest possible remnant fields are given). As in 
the RSO mode of measurements sample travels only in a 
very central part of the magnet bore, we can treat this 
field as an additional homogeneous field acting on the 
sample, but with very important characteristics: it is di- 
rected opposite to the last experienced strong field by 
the magnet and assumes the largest values around zero 
magnetic field [32]. This results in a well known "neg- 
ative" hysteresis, an apparent inverted (clockwise) mag- 
netic field loop seen in soft ferromagnets when the results 
are plotted against Hj, the field calculated from the cur- 
rent sent to the magnet and reported in the MPMS sys- 
tem's data files. This fictitious effect occurs only when 
this nominal Hj is used instead of the real field Hr taken 
from a field sensor located in the close vicinity of the sam- 
ple (which does not exist in the MPMS). Since in this 
chapter we consider the effects caused by the difference 
between Hj and Hr, we will strictly stick to this nota- 
tion here. We abandon this distinction in the rest of the 
text and use the symbol H which actually indicates Hj 
(values given in the .dat files) unless it is clearly stated 
otherwise. 

The very existence of H tr (resulting in a discrepancy 
between Hj and H^) and the use of Hj results in an ap- 
parent ferromagnetic-like response seen in samples which 
the magnetic response is dominated by a diamagnetic 
contribution. As the example of such an effect we dis- 
cuss a room temperature hysteresis loop m(Hj) collected 
for a sapphire substrate. It is important to note that 
there may exist a 'genuine' FM component in sapphire 
14] and in particular in a non-properly handled sapphire 
27| . We discuss this issue further in this review (chapter 
IIII ). Here we limit ourselves to magnetically clean mate- 
rial, which has been handled in the best possible way to 
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FIG. 6. (Color online) Left axes: nonlinear component with 
respect to Hj, the magnetic field reported in MPMS *.dat files 
which is calculated upon the current sent to the magnet, as a 
function of Hj for (a) 60 mg diamagnetic sapphire substrate 
and (b) paramagnetic salt Dy2 03. The degree of compensa- 
tion of original m(Hj) exceeds 99.9% above 10 kOe in both 
cases. Please note the opposite chirality of the two reported 
Am(Hj). The right axis in (b) denotes calculated values of 
the trapped field in the magnet, see the text. The dashed 
grey lines indicate the switching field when the power supply 
is switching from high resolution to low resolution charging 
mode. 

avoid any magnetic contamination. The actual specimen 
for this test consists of two 5 x 5 x 0.3 mm 2 pieces of total 
mass of 60 mg. The direct measurement gives fairly linear 
in Hj values of the order of 2 x 10 -5 emu/kOe, so in or- 
der to check whether a ferromagnetic-like contribution is 
present we compensate the MPMS data by a linear, aHj 
term, with a manually adjusted to flatten out the result- 
ing dependency in a widest possible range of Hj. There- 
fore a can be regarded as a specific susceptibility of this 
specimen, xs- The resulting Am(Hj) = m(Hj) — aHj is 
presented in Fig. H{a). Although resembling a ferromag- 
netic response, such a shape is caused by the presence 
of the trapped flux H tr in the superconducting magnet. 
Other evidence for the dependence of this type of resid- 
ual hysteresis in net-diamagnetic samples on the history 
of the magnet has been provided before [13] . Here it will 
be demonstrated that identical residual hysteresis with 
opposite chirality can be recorded when paramagnetic 
samples are used. 

To substantiate the origin of the residual hysteresis in 



sapphire and to determine the size of the residual field of 
the superconducting magnet we turn to the Dy2C>3 sam- 
ple (a paramagnetic salt) provided by Quantum Design 
for every MPMS system equipped with AC option[3J]. 
We measure this test sample with the identical experi- 
mental procedure (i.e., identical history of the magnet) 
as the sapphire substrate and compensate for the lin- 
ear term, f3Hj, this time adjusting j3 to obtain Am(Hj) 
symmetrical with respect to Y axis at high fields , see 
Fig.lHlJb) dark blue points. Importantly, this time we note 
the clockwise character of Am(Hj) loop (an inverted hys- 
teresis), the opposite one to that of the sapphire. Using 
the well-established value of (3 we can calculate the real 
magnetic field acting upon the sample Hr = m{Hj)/f3 
and obtain the "field error" H Tr = AH = H R -Hj. The 
values AH are indicated on the right-axis in Fig. [BJb). 
To cross-check the validity of our approach, we repeat 
the same procedure to calculate AH for the sapphire 
substrate taking this time slightly altered value a as xs. 
We add the results to Fig. [SJb) (marked in grey) and the 
full equivalence of both sets of AH confirms the common 
origin of the observed Am: the extra field acting on the 
specimen, the H tr field. The different chirality of both 
Am{H) loops can now be readily understood, as H tr is 
the same positive or negative field for the sweeps up or 
down, respectively, acting on substances with an opposite 
sign of the susceptibility. 

On the other hand it is very instructive to note that we 
needed a marginal change of a (less than 0.05%), from 
2.1465 x 10~ 8 emu/Oe used to calculate Am in Fig. [6](a) 
to 2.1457 x 10~ 8 emu/Oe used to asses AH in Fig. Mb) 
to effectively remove the deceiving ferromagnetic-like ap- 
pearance of Am presented in Fig. Eta). It shows how 
easy it is to falsely interpret the data. It is therefore in- 
structive to see what happens when one disregards H tr 
and considers such Am from Fig. [SJa) as a genuine fer- 
romagnetic response. In this case the 'saturation' mj 
(or a 'remanence', if one considers the H — case) is 
m ~ 3 x 10~ 7 emu. Let's imagine now that this is not 
a pure test-substrate but a typical sample consisting of 
a thin paramagnetic DMS layer deposited on a similar 
sapphire substrate. A recent investigation of a; < 1% 
(Ga,Mn)N is a good example of such a study, see [20j. 
Taking a typical thickness of the DMS layer of ~ 200 nm 
one easily notices that S, the value of spin per each Mn 
is proportional to l/x(%). Therefore, calculating S per 
Mn ion from the artificial of sapphire, S exceeds the 
maximum expected value for Mn (S = 5/2) already at 
room temperature in such layers with x < 0.4% even if 
the Mn is completely unpolarized. Remarkably, the effect 
will be stronger for thicker substrates, thinner layers and 
in nominally stronger magnets {H tr approaches 30 Oe for 
the 70 kOe MPMS). 

On the contrary, a presence of a strong paramagnetic 
component in the whole signal of the investigated speci- 
men, frequently met in DMS layers at low temperatures, 
will reduce either absolute values of the overall m at low 
magnetic fields or a coercive force (if present). It also 
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leads to 'negative' hysteresis loops if the true coercive 
force is in the few-Oersted range. More interestingly, the 
-ff tr -induced effects will be strongly temperature depen- 
dent, and so they may create a number unusual and 'in- 
teresting' findings. In any case, depending on the history 
of the field in the magnet, an iJj-based low field evalua- 
tion of any magnetic signal will be at least inaccurate, if 
not completely wrong. 

Since the flux pinning is a property related to the type 
II superconductivity of the NbTi, the basic material from 
which the superconducting wire of the magnet is made of, 
the only way of minimizing H tr is by warming the magnet 
above 9 K, the superconducting critical temperature of 
NbTi. As in general this is highly impractical idea, there 
is a very useful option available, called "magnet reset", 
which, upon our experience in weak signal magnetism, we 
regard as a kind of "must have" for everyone interested 
in sensitive and precise magnetometry, at weak magnetic 
fields in particular. This option features one layer of 
heater wire wound along the entire length of the magnet 
and radially about half way through the magnet. A heat 
pulse through this heater will locally turn some windings 
into the normal (resistive) state which will then dissipate 
more energy (heat) to drive some neighboring windings 
to the normal state, and so on. Finally the whole super- 
conducting material is driven into the normal state, and 
so the all pinned flux is expelled from the body of the 
magnet on an expense of a little amount of liquid helium 
in otherwise helium full dewar. It should be noted that 
this option is only useful in cases when a small magnetic 
field is still in the magnet. At nominally zero magnetic 
field the use of the reset option may not heat the entire 
magnet and — opposite to what is intended — H tr may in- 
crease or even reverse sign, see [13] > Fig- 6. 

Unfortunately, the magnet reset can purify the mag- 
net exclusively for operation around H = 0. It is of 
no help during typical m(H) measurements, particularly 
those requiring the full extent of the magnetic fields. In 
such a case, we see no other method to eradicate H tr 
from post measurements analysis then by numerical re- 
placement Hj in data sets by Hr using H tr . To make 
it possible one has to establish the Htt(Hj) dependence 
for the specific system and the field (and eventually tem- 
perature) history of the magnet. 

The bottom-line is that there is no universal 'table' of 
Hxr(H) as H tr depends on the most recent magnet his- 
tory. This is illustrated in Fig. [JJ We depict there H tr in 
±10 kOe window established during a multiple field loop- 
ing up to ±50 kOe through progressively stronger mag- 
netic fields in the MPMS XL 5 T in Warsaw, Fig. [TJa) 
(data presented in Figs. [BJb) and[9ja) are a part of this 
large data set). We started having the magnet purified 
with a self written degaussing sequence followed by the 
magnet reset. What can be immediately seen is a change 
of the initial 'slope' of the data in the field range to 
±1 kOe (dark blue) compared to the rest of the field 
sweeps. Importantly the slope is preserved as the field 
is limited to this range and then it starts to gradually 
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FIG. 7. (Color online) (a) A development of the remnant 
(trapped) field in an MPMS 5T magnet studied through 
sequential field sweeping to progressively stronger magnetic 
fields showed in ±10 kOe window. The initial part, ±1 kOe 
and to 2 kOe (dark blue) was measured after nulling the field 
by using with a self written degaussing sequence followed by 
the magnet reset. Then the field was swept to -3 kOe (light 
blue), 5 kOe (dark green), -5 kOe (green), 10 kOe (red), - 
20 kOe (brown). The remaining colors indicate field sweeps 
from: -20 kOe (orange), -20 kOe (orange), 32 kOe (yellow), 
-32 kOe (magenta), 50 kOe (purple), and -50 kOe (grey). The 
remnant field was calculated upon magnetic moment data of 
paramagnetic Dy2C>3 described in the text. Arrows of the 
corresponding colors indicate the field sweep directions. The 
dashed grey lines indicate the switching field when the power 
supply is moving from high resolution to low resolution charg- 
ing mode, (b) Dark blue: the lowest field part of the data from 
(a). The grey marked points were collected after the whole 
sequence from (a) followed by magnet reset option. Note a 
small negative field, Hr ~ —0.2 Oe, still present in the mag- 
net. The data for this panel were compensated using slightly 
different factor to level them off at ±1 kOe range. 



change during the following field sweeps to 2, -3, 5, 
and so to ±50 kOe. Importantly, if the same field loop 
is performed twice (not shown here), the obtained H tr 
values are very much the same. So, if a particular type 
of measurement is performed frequently, its pays off to 
established the specific set of HT r {Hj) and use it to cor- 
rect Hj with a great confidence. It is worth noting that 
this method also automatically compensates for a slight 
change of Hj at ±5 kOe. A small jump in the data is 
seen in Figs. [5] and [7] and the relevant switching fields are 
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marked by vertical grey dashed line. These jumps occur 
when the power supply is switching from high resolution 
to low resolution charging mode, see chapter fllC 11 

A necessity of performing the degaussing procedure to 
establish a true zero field in the magnet is exemplified in 
Fig-E^b) where grey points mark H tr measured after per- 
forming the magnet reset alone after the described previ- 
ously field sweeps to 50 kOe. For this panel we compen- 
sate the data by slightly different value of /3 (increased by 
0.6%) to blow up a change of Ht r during field sweeping in 
weak field range. We note that the 'grey' Ht r (Hj) values 
are shifted down with respect to the initial ones, and the 
overall spread of the points (marginally weak hysteresis) 
covers the range of less than 0.5 Oe at fields below 1 kOe. 
This indicates that the absolute field reproducibility at 
these weak magnetic fields can be maintained better than 
0.05%. 



3. Flux creep and escape 



1 1 1 

■ a) after 10 
L no over 

1% 


i ■ i I i ■ i I ■ ■ ■ I ■ i ■ 

kOe GaAs T = 300 K 

shoot 

^ H = 1 kOe 

b) after 10 kOe c) after 40 kOe 
and and 
magnet reset magnet reset 











2 4 6 8 10 

time ( h ) 



FIG. 8. (Color online) Field drift in 5 T magnet of MPMS 
system. The magnetic response at Hj = 1 kOe was recorded 
each time for approx. 3 hours for a piece of GaAs substrate 
after: (a) the field was brought down directly (in the no over- 
shoot mode) from 10 kOe, (b) after the initial 10 kOe the 
magnet was firstly purified using the magnet reset option and 
then the 1 kOe was set, and (c) the same as (b), but with 
the starting field of 40 kOe. The negative H tr of the order 
of -15 Oe present in (a) is responsible for the vertical offset 
between (a) and (b-c). 



Even in the persistent mode of the magnet, some of 
the pinned magnetic flux lines in the superconducting 
wire will redistribute with time and some will leave the 
wire material. These that enter the bore of the magnet 
must be represented by an extra nonvanishing current in 
the magnet, so the magnetic field will not be stable due 
to the very nature of the type II superconducting mate- 
rial of the magnet (for an extensive description see (32|). 
This property has got a profound implication in low field 
magnetometry, e.g. when one wants to establish x(T) of 
a paramagnetic nanostructure at high temperatures (see 
[35j). Therefore a subtraction of m(T) of the substrate 
measured at the same H (normalized by the correspond- 
ing sample weight) must be made with the best possible 
accuracy. The time instability and history dependence 
of the magnetic field in the magnet renders this task im- 
possible if one does not take the flux creep into account. 
The exact value of Hr is history dependent and if one 
sets a weak filed H = Hj the obtained Hr may be up 
to ±15 Oe off the target (see Fig.JTJa)). One can try to 
live with this by sticking to exactly the same experimen- 
tal routine, assuring repeatable Hr, but the flux creep 
and escape cause the whole m(T) to be performed at 
a continuously varying magnetic field, as is indicated in 
Fig. [H part (a). There the field was brought down to 
1 kOe from 10 kOe in no overshoot mode resulting in a 
field drift on the magnet over hours. In this particular 
case, the magnet reset option was of great help (as indi- 
cated in parts (b) and (c) of the same measurement) by 
purging the magnet from the trapped flux directly be- 
fore switching the 1 kOe field on. Thus it allows not only 
for setting a perfectly repeatable field values but also it 
sets conditions for an order of magnitude smaller field 
drift, independently of the most recent field history in 
the magnet. 



4- Limits of high field resolution 

As already pointed out, accurate studies in nanomag- 
netism require a great deal of compensation to be in- 
volved and obviously, the accuracy of a difference of two 
comparable numbers directly depends on the absolute 
precision of both the minuend and the subtrahend. As 
the system establishes m on the account of a numeri- 
cal comparison of V(z) to T(z) [36| then any deviation 
of the functional form of V(z) from that of T(z) will 
yield somewhat distorted value of m. We routinely de- 
tect such distortions in all the scans obtained at magnetic 
field stronger than about 20 kOe. 

To visualize the effect we present in Fig. Eta) compen- 
sated for a linear in magnetic field component room tem- 
perature data of the same piece of the paramagnetic salt 
(Dy 2 Os) as discussed previously in section Hi C 21 Here 
we focus on the sizable outturn of the Am(if) from the 
expected flat dependency on H seen at above 25 kOe. 
This nonlinear component in the original m{H) cannot 
be simply related to deviations from an ideal linearity of 
this particular paramagnet, as qualitatively the same ef- 
fect is seen for diamagnetic specimens, as exemplified in 
Fig. E^b) . Neither this effect can originate simply from 
an excess of noise at high fields, since it is seen for both 
strong and rather weak signals and is smooth and mono- 
tonic on H. Nor it is caused by the presence of a remnant 
field in the system, as contrary to Ht r (H), the amplitude 
of the deviations grows with the field and it is at least an 
order of magnitude stronger than reduction of the mag- 
netic signal due to the presence of H tr - Then, the last 
obvious cause would be a systematically growing uncer- 
tainty of the larger magnetic fields in the system. In 
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FIG. 9. (Color online) (a) and (b): (brown bullets) nonlin- 
ear in H component of m for Dy2C>3 paramagnetic salt and 
sapphire, respectively. The arrows indicate the direction of 
field changes during these measurements. The compensation 
method and coefficients a and ft are the same as described 
in chapter III C 2 1 (c-f) (dark green bullets) distortions scans 
ri(z) at 5, 20, 40, and 50 kOe with respect to a scan collected 
at 10 kOe (scaled to the relevant field). Please note that the 
vertical scale for T){z) is in millivolts, not in volts as it is typi- 
cally for full V(z) indicating an enormous compensation level 
of V(z). Red lines: T(z) fitted to rj(z) in the linear regression 
mode. Results of relevant fittings for all the experimental 
points are marked in red in (a-b). They cross Y — line at 
H = 10 kOe as this is the measurement that all the points 
were compensated against. The green lines there are the same 
as red, but with slightly adjusted slope to level them off at 
moderate fields. They match then the original Am perfectly, 
underlying a remarkable linearity of the regression routine. 



order to check this possibility we numerically compen- 
sate the raw V(z) scans to asses the symmetry of the 
resulting distortion scan rj{z). To do so, we take one of 
the scans from the undistorted part of the m(H) as our 
reference scan (V(z) obtained at H = 10 kOe was se- 
lected for this purpose). Then we scale this V(z) by the 
relative value of the magnetic field of each experimental 
point and subtract such a scaled reference accordingly to 



obtain sought distortions r](z) [33]. Four examples of r)(z) 
are established this way and are presented in Fig. c-f). 
They clearly indicate an odd character of the distortions 
at large fields and so decisively ruling out any inaccura- 
cies related to a systematic difference between the real 
magnetic field acting on the sample and the field values 
reported in the data files. 

At the moment we cannot identify the effect responsi- 
ble for the observed deviations at higher magnetic fields. 
We have already pointed out (see chapter III B 3[) that an 
odd component in V(z) is a fingerprint of a presence of a 
somewhat tilted magnetic moment. Although very likely 
at weak magnetic fields, here, the moment-tilting should 
be ruled out based on the increases of the directional uni- 
formity of the field generated in typical solenoids with 
increasing current in the magnet. Even if the samples 
would possess any magnetic anisotropy, it should nei- 
ther be visible at strong fields nor increase with the field. 
However, whatever the underlying reason, the presence 
of this odd r](z) that causes the least-squares regression 
routine to be best satisfied for slightly lower values of m 
than it would be for a 'clean', a perfectly Y(z)-like V{z). 
This can be confirmed by fitting T(z) to the distortion 
i](z). The fitted T(z) are marked in panels (c-f) and the 
corresponding moments (for all the experimental points) 
are plotted in Figs. |Ha-b) in red. Their different slope is 
due to differently levered compensation (giving Arii = 
at H = 10 kOe this time), but after re-adjusting the 
slope the moments obtained from fitting rj(z ) m atch the 
original results nearly perfectly (green lines) [38| . As the 
numbers in Fig. [9] indicate, this 'property' effectively pre- 
cludes reliable determination of the third significant digit 
of the magnetic moment by the MPMS at fields stronger 
than 25 kOe, and so the compensation methods are lim- 
ited to about 99-99.5% of compensation of the total signal 
at high magnetic fields. Finally we observe that despite 
the fact that the compensation exceeds 99.5%, the perfect 
match between values obtained cither by a compensation 
of the results of the fitting (the MPMS results) or by a 
fit to compensated raw scans (both compensation by a 
linear term in H) indicates a truly remarkable linearity 
of the regression routine. 



5. Auto-range option 

It has been reported before that the auto-range option 
of the SQUID magnetometer can lead to artifacts once 
the emu-range is changed upon increasing or decreasing 
total magnetization signals and the characteristic signa- 
tures of this type of artifacts have been described in de- 
tail before [231. Here it should be noted that the emu- 
range-artifact is not present any more is that particular 
machine (the MPMS in Duisburg) whereas in the other 
system in question in the present work, it has never been 
observed. On the other hand, it was observable in data 
recorded with other machines as well, e.g. in Ref. (l7j . 
Fig. 1. The disappearance of the emu-range artifact in 
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the MPMS was correlated with the following three ac- 
tions which unfortunately were carried out at the same 
time: exchanging of the yellow cables connecting the rf- 
SQUID-amplificr with the cryostat, re-tuning the SQUID 
and re-adjusting the filter settings of the rf-amplifier. In 
any case, this indicates that this particular artifact is 
indeed correlated with an improper capacitance-match 
of the R-C-circuits responsible for the various " range" - 
values as discussed before [27|. Once such artifacts are 
noticed, it is therefore advisable to re-tune the SQUID 
and to check the filter settings. In the worst case, it 
should be considered to replace the yellow cables. At 
this point it should be noted that the feed throughs at 
the cryostat are in some cases not He-leak tight and He 
blow-off may result in cold plugs and thus condensation of 
water in these connectors which in turn alters the capaci- 
tance of the cables. An " 'easy fix" ' is to fill the connectors 
with vacuum grease. Unplugging the yellow cables and 
drying them out with a heat gun is less advisable, since 
the connectors are specified for only 40 cycles of plugging 
and unplugging. 



D. Sample mounting 

The (second order) gradiometer configuration of the 
SQUID detection coils, that is the equal number of turns 
in each directions arranged as shown in Fig. [2] has got 
two profound implications on the magnetomctry. Firstly, 
the pick-up assembly is far less sensitive to distant mo- 
ments than those within the volume of the gradiome- 
ter. Secondly, there would be no magnetic signal from 
an infinitely (= sufficiently compared to the size of the 
gradiometer) long object, providing it is magnetically 
uniform, independently whether it moves or not and 
what kind of magnetic response dominates its proper- 
ties. Practically, anything of a length of 15-20 cm what 
allows to fix the sample close to its center and has got an 
adaptor connecting it with the main (long) sample rod 
will work for sensitive magnetometry, however, the least 
'magnetic', the better. The MPMS XL system manufac- 
turer advices to use clear drinking straws as sample hold- 
ers. Indeed, they are of generally low magnetic signature 
(the background-signals usually stay below 10 -7 emu) 
and provide a great versatility when sample fixing is con- 
cerned. Since typical specimens are usually of a size of 
3^5 x 5.5 mm 2 , all can be held in place inside the 
straw without any other means than by clamping them 
in-between the two walls of the straw as illustrated by 
Fig. HUk). In this configuration the background of the 
sample holder is minimal and the magnetic field is ap- 
plied in the sample plane. To reliably perform in-plane 
versus out-of-plane measurements it is necessary to have 
a sample piece which can be freely rotated, e.g., by us- 
ing two wooden sticks as non-magnetic " manipulators" , 
inside the same straw. For that purpose, a~4x4 mm 2 
sample piece can be clamped into the straw along its di- 
agonal as illustrated by Fig. HUb ) . Both options assure 



sufficiently well-centered samples with respect to the gra- 
diometer axis if straight straws are used. 




FIG. 10. (Color online) Different means of sample mounting 
for SQUID measurements using a straw: a) in-plane mea- 
surements for sample sizes 3 — 5 x 5.5 mm 2 , b) in-plane ver- 
sus out-of-plane measurements allowing to rotate the sample 
(~ 4 x 4 mm 2 ), and c) solution for smaller sample sizes using 
a second straw. 

In cases where the sample piece is to small to be held 
freely inside the straw (less than 4x4 mm 2 ), one has to 
use a second straw which has to be carefully slitted along 
its length ideally using a clean ceramic blade. This sec- 
ond straw can be slightly twisted up while inserted into 
the outer straw and presses the sample piece to the wall 
of the outer straw as illustrated by Fig. [TOb). The last 
configuration has to be used in a greater caution when ac- 
curate absolute values are desired, as the substantial off- 
centering of the whole specimen will increase the strength 
of the coupling with coils (see chapter lllB 2j) . In any case, 
one has to strictly avoid cutting any holes into the walls of 
the straw, or even severely denting/puncturing the straw 
in the region of the sample. This inevitably creates an 
additional (temperature-independent) signal which can 
be effectively diamagnetic (accumulation of material) or 
paramagnetic (lack of material). 

Having all the above-mentioned precautions in mind 
the straws may be of a great value for everyday usage, 
even precise, magnetometry in a range down to a few 
10 -7 emu. For such purposes and to reduce the danger 
of external magnetic contamination one should avoid re- 
using straws as well as refrain from cleaning the straws 
prior to the measurement. In cases where really small 
signals (in the 10 -8 emu range) are to be researched, one 
has to look for an alternative. In fact, about 50 straws 
from two different batches provided by Quantum Design 
were tested. They were scanned along almost their entire 
length [39] at 3, 1 kOe and at remanence at 10 and oc- 
casionally at 200 K. Then the collected long scans were 
recalculated to obtain a corresponding (fictitious) mo- 
ment distribution along their length, the magnetic profile 
fi(z). All the test were done without 'front' (top) and the 
end straw adapters fitted. We wanted to test the straws 
themselves. The presence of straw adapters produces so 
called "paramagnetic Curie tail" at the lowest temper- 
atures and relatively strong fields [loj]. We found that 
the typical amplitude of /d(z) was a few times 10 -7 emu 
at 10 K and 1000 Oe, and only in two of all the tested 
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FIG. 11. (Color online) Magnetic moment distribution fi(z) 
at the central part of the three different sample holders. The 
grey line represents one of the very best straw from those 
provided by Quantum Design, an average straw has got 2- 
3 times larger amplitude of n(z). The dashed orange line 
represents a typical response of a 0.5 mm gold wire, while 
the black line gives a typical fi(z) distribution in 2 mm wide 
silicon strip. The position z is counted from one of the ends 
of these sample holders. 



straws fi(z) was below 10 -7 emu, see an example given 
in Fig. 111! Symptomatically, we found straws from one 
'box' to be substantially " less magnetic" (i.e. more mag- 
netically uniform) than from the other. Finally, some 
proved to be very heavily contaminated, despite the fact 
that all the straws were thoroughly cleaned in organic 
solvents prior to their tests. So, to cut on this tedious, 
time and resources consuming testing job, the Warsaw 
team started to use a similar length of a straight 4N sil- 
ver or gold 0.5-0.7 mm wire. They exhibit as a rule the 
same small /i(z) as the best straws found by us (dashed 
orange line in Fig. 1111) . are easy to maintain in 'pristine' 
condition for a very long time. There are, however, two 
disadvantages related to this kind of a sample holders. 
Obviously, a glue is needed to fixed the sample, secondly, 
the wires deform or bent easily. The last issue has been 
solved by switching to 1 to 2 mm wide silicon strips pro- 
fessionally cut from 8" commercial wafers 41]. They are 
fixed between the standard straw adapters and so are 
easy to assemble and provide mechanical stability dur- 
ing the measurements. Our tests show that all silicon- 
based sample holders conform to our stringent require- 
ments concerning magnetic uniformity as we find their 
n{z) to be practically equal to the noise level of the mea- 
surement, see black line in Fig. QTJ However, the major 
drawback is again the necessity to using a glue to fix the 
sample. After testing many substances for a possible con- 
tribution to the measurement as well as for mechanical 
reliability and low temperature behavior we settled for 
the new branch of the GE varnish [42[ (heavily diluted), 
which proved to best fulfill all requirements. It must be 
finally noted that the usage of Si strips or similar bodies 



may lead to a sample off centering, so it should be taken 
into account in very precise magnetomctry, see chapter 
lHB2l 



III. ARTIFACTS DUE TO IMPROPER SAMPLE 
HANDLING 

It is needless to say that every object of any investiga- 
tions should be as well defined as possible and free from 
interfering alien contaminants. This is particularly true 
in integral (SQUID) magnetometry where the typical sys- 
tems integrate the magnetic response from the specimen 
and its vicinity. The issue has become of a prime im- 
portance how the vast part of the community switched 
to magnetic investigations of nanoscale magnetic objects 
and this is there where the main advantage of SQUID 
magnetometers, their enormous sensitivity, has become 
also the main concern. This is clearly seen when one re- 
alizes that a single 1 ng weighting 5 ^.m swarf of iron will 
produce a moment of 2 • 10~ 7 emu, and a typical flake of 
rust floating in the air around can be easily a few times 
larger. Thus, the presence of a virtually undetectable 
(even using a decent optical microscope) contaminations 
can totally overwhelm the actual response of a small- 
signal sample. Some contaminations produce a response 
which is indistinguishable from that of the sample, so a 
great care must be exerted at every experimental step, 
starting from the choice of a substrate, through the de- 
position and processing until the final mounting on the 
sample holder. But this is usually not enough, so a very 
extensive testing augmented by some experience in the 
field is required to trustworthy separate the real signal 
of interest (if any) from spurious effects. However, occa- 
sionally, contaminations may reveal themselves through 
a unique, even bizarre, response, effectively very soundly 
alarming the experienced investigator. We believe that it 
is very instructive to consider one of such cases in greater 
detail. 

Such a "fortunate" case is exemplified in Fig. [T^] which 
presents results of the most obvious and presumably com- 
mon mistake made during the preparation of a sample 
for SQUID measurements, namely an edge contamina- 
tion by a (ferrous) material stemming from tweezers or 
another tool used to prepare, cut, or handle the specimen 
[T3 . [l3j . This issue has become of a great importance 
when many laboratories started to investigate materi- 
als deposited on substrates made from some of hardest 
substances like sapphire or SiC which are very likely to 
abrade the softer material of the tools they are handled 
by. But even Si requires a great deal of care because of 
its strong tendency to bind iron or nickel [lH, |43| , or to 
absorb on its surface submicron Fe-rich particles during 
etching from even a high grade KOH j44j. Here we con- 
sider the case of a sapphire substrate and in Fig. [T2Ta) 
an m(H) measurement at 300 K of an as-cleaved GaN 
layer on a sapphire(OOOl) substrate is shown after the 
diamagnetic background has been subtracted. A clear 
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FIG. 12. (Color online) Hysteresis loop of an as-grown GaN 
layer on sapphire(OOOl) recorded at 300 K in the as-cleaved 
state (full symbols) and after cleaning in HC1 (open sym- 
bols), (a) shows the data after the diamagnetism has been 
subtracted, (b) the low field part of the full signal, that is 
with the diamagnetic signal of the substrate, and (c) shows 
the respective fitted z-position of the SQUID signal (iterative 
regression mode). 



ferromagnetic-like hysteresis is visible (full squares) with 
one unusual but alerting feature, namely a kink seen be- 
tween 10 mT and 25 mT. The first important hint to 
unravel its origin comes from the inspection of the full 
moment of the sample, the uncompensated m(H), pre- 
sented in Fig. IT2T b). We notice that the kink is located 
in the field range where the raw data changes sign, that 
is around the compensation field He for which values 
the diamagnetic background become comparable to that 
of the "ferromagnetic" signal. The second and the most 
important clue comes from Fig. IT27 c) which shows the 
fitted "best" z-positions of the individual SQUID mea- 
surements obtained in the iterative regression mode of 
the MPMS. It is obvious that the kink corresponds to 
data points where the fit returns an apparent sample po- 
sition which is up to 7 mm shifted with respect to the 
actual position. This can be understood by considering 
the following: the ferromagnetic signal originates from 
the edge of the sample at the top (or bottom), whereas 
the diamagnetic signal comes from the entire sample, i. e., 
it stems from the center of mass of the sample. This ex- 



perimental situation is depicted in Fig. [3jb). At mag- 
netic fields where both signals are of approx. the same 
size, the respective curves recorded by the MPMS have 
an JV-shaped appearance rather than the expected field 
distribution with a clear extremum in the center. Thus, 
the fitting routine picks one of the two extrema of the JV- 
shaped curve and returns a position which is off-center. 
It is characteristic that this off-set changes sign because 
the fitting routine takes the extremum of the JV-shaped 
curve which is initially bigger. If the other extremum of 
the iV-shaped curve becomes larger at a different field, 
the fitting jumps to the other extremum of different sign, 
therefore causing the characteristic change of sign in the 
z-position plot as shown in Fig. [T2T c). Characteristically, 
at strong fields the z-position goes to zero, as the prop- 
erly centered diamagnetic signal dominates, whereas at 
H = z ~ 3 mm, thus very close to the edge of the 
sample where the contaminant was. And in both lim- 
its the shape of the measured V(z) should be very much 
like the ideal one, which contrasts the case of the tilted 
magnetic moments (see chapter III B 3|) where neither the 
iterative nor the linear mode of fitting will be acceptable. 
One may trace the regression fit value for this (reported 
by the system in .dat files), we would rather advice to 
inspect the raw data to visually asses each 'problematic' 
case. 

The occurrence of a little kink in the hysteresis at zero 
measured total magnetization together with the charac- 
teristic up-down behavior of the z-position in the itera- 
tive regression mode as shown in Fig. [12] is thus a clear 
indication for ferromagnetic contaminants at the edge of 
the sample. Note that these features would be not visible 
in the data, if the magnetic contamination is not located 
at the top or bottom but at the left or right edge (close 
to the center) with respect to the movement direction of 
the sample. The safest procedure is to measure m(H) 
curves with two by 90° rotated in-plane configurations. 
A comparable artefact occurs if a non-homogeneous fer- 
romagnetic sample on a diamagnetic substrate is inves- 
tigated, e.g., nanoparticles deposited at one side of the 
substrate, ion implantation of only half of the sample or 
lithographically prepared structures with uneven distri- 
bution. Furthermore, if one uses the linear regression 
mode of the MPMS - where only the amplitude but not 
the position is fitted - one would lose any experimental 
evidence for magnetic contaminations from the edges of 
the sample and the data in Fig. IT2Ta) would roughly fol- 
low the red line. That the apparent ferromagnetic-like 
signature in Fig. [12] is indeed caused by a ferromagnetic 
contamination of the edges of the sample is confirmed by 
etching the sample in HC1 followed by thoroughly clean- 
ing the sample in acetone, ethanol and ultrapure water 
in an ultrasonic bath for 5 min. each. The M(H) curves 
were measured once again (open circles) and it is ob- 
vious that the magnetic hysteresis has almost vanished 
and the kink has disappeared (the marginal opening orig- 
inates from not compensating the data for the presence of 
the trapped field in the magnet, see chapter Til C 21) . The 
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up-down behavior of the fitted position is now found at 
zero magnetic field, where the total signal is zero and it 
is much less pronounced. The overall magnetic behavior 
is now purely diamagnetic. Furthermore, the role of the 
GaN layer can be assessed by comparing these findings to 
the ones in Ref. [13], where a similar behavior has been 
reported for bare sapphire. 

The magnetic contamination at the edge of the sample 
most likely originates from cleaving it into suitable pieces. 
It is very common to use a diamond stencil for that pur- 
pose; however, the diamond is usually a tiny piece at- 
tached to a stainless-steel rod. If this rod touches the 
edge of a very hard substrate material it is very likely 
that small particles of the metal are abraded. Hence a 
thorough cleaning of the sample, especially if sapphire 
or SiC is used is crucial. In case of substrates which are 
much softer than stainless steel, the situation is less criti- 
cal and cleaving usually does not produce any detectable 
ferromagnetic contamination. Of course no stainless steel 
tweezers should be used, see Ref. [l2| in any case. Teflon, 
Titanium or Zirconia tweezers are suitable but one has 
to assure that they are not used, e.g., for nanoparticle 
samples, in-between since the tweezers-handling always 
transfers small amounts of nanoparticles to the edges of 
the sample. 

Once great care has been taken to avoid magnetic con- 
tamination of the sample, a sequence can be performed 
to asses the intrinsic performance of the SQUID for an 
"every-day" measurement. For that purpose and to "sim- 
ulate" straightforcward standard measurements without 
tedious precautions and correction procedures, the iden- 
tical sample (including straw) as in Fig. [12] has been 
further measured at low temperatures following a stan- 
dard sequence of recording an m(_ff)-loop at 300 K, then 
at 5 K and then two m(T) measurements while warm- 
ing under field cooled (FC) and subsequently after zero- 
field cooled (ZFC) conditions. The outcome is plotted 
in Fig. [T3] after subtracting the diamagnetic background 
which was determined from the mid-field (6 to 10 kOc) 
m(77)-behavior at 300 K, the most simple background 
correction possible. Figure ITBT a) shows two hysteresis 
curves measured at 300 K and 5 K, revealing virtually 
identical magnetic response at low temperatures and high 
magnetic fields with a '"saturation"' magnetization of 
4 • 10~ 7 emu. At low fields the residual hysteresis be- 
comes visible at 300 K (b). A "remanence" of 2-10 -7 emu 
and a rounded loop with "coercivity" is visible as well. 
At 5 K roughly the same "remanence" and "coercivity" 
is visible. These small residual signals are very simi- 
lar to the findings in Fig. [HJa). Since the GaN/sapphire 
sample is net-diamagnetic as well, it is highly likely that 
the findings are dominated by H tr in the superconduct- 
ing magnet. The absence of a true ferromagnetic sig- 
nal for this sample can be independently inferred by the 
FC/ZFC curves in Fig. ri"3T c) which are not separated at 
any temperature within the noise level and therefore are 
indicative of the absence of any ferromagnetic or blocked 
superparamagentic response of the sample. Note, that 
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FIG. 13. (Color online) (a) Residual magnetic hysteresis of 
a GaN layer on sapphire(OOOl) substrate at 300 K (squares) 
and 5 K (circles). No additional paramagnetic signal can be 
seen at low temperatures and high magnetic fields, (b) The 
apparent hysteresis is not temperature dependent and field 
cooled (squares) versus zero field cooled (circles) curves coin- 
cide indicating the absence long range ferromagnetic order or 
blocked superparamagnetic behavior (c). 



the m(T) measurements were performed while sweeping 
(no stabilization of the temperature for each data point) 
T at a high rate (7 K/min) with little statistics for each 
point (1 scan, 3 cycles). Even under these "everyday" 
conditions the residual artifacts as discussed here never 
exceed 4 • 10~ 7 emu and relative changes of the order of 
1 • 10~ 7 emu can be resolved. However, this is only true, 
if proper centering according to Fig. [TUT a) is used and 
magnetic contamination is ruled out. For users with lit- 
tle access-time to SQUID magnetometry this may be the 
ultimate performance achievable and reliable detection of 
signals below 4 • 10 -7 emu require much more time, care 
and efforts. Nevertheless, this "realistic" sensitivity still 
makes the MPMS one of the most sensitive commercial 
magnetometers . 



IV. CONCLUSION 

In conclusion our detailed analysis of artifacts and 
limitations of commercial MPMS SQUID magnetome- 
ters demonstrate that reaching the intrinsic sensitivity of 
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about 5 TO -9 emu (according to the manufacturer's spec- 
ification) is very challenging if not impossible for samples 
of finite size, where various magnetic constituents sum up 
to an overall magnetic response of the specimen. Even 
for magnetically homogeneous samples, reaching the ul- 
timate sensitivity is difficult when tilted moments occur 
or the sample is slightly off-axis. In addition, the are 
artifacts which originate from the intrinsic principle-of- 
operation of this gradiometric setup. Finally, the use 
of a superconducting coil to generate the magnetic field 
adds more complexity, especially for magnetometry at 
low magnetic fields. Where possible, we have tried to 
provide typical signatures, concepts to circumvent prob- 



lems or to compensate for artifacts. 
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